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The three-band Emery model is reduced to a single-particle quantum model of Falicov-Kimball 
type, by allowing only up-spins to hop, and forbidding double occupation by projection. It is used 
to study the effects of geometric obstruction on mobile fermions in thermodynamic equilibrium. 
For low hopping overlap, there appears a plateau in the entropy, due to charge correlations, and 
related to real-space disorder. For large overlap, the equilibrium thermopower susceptibility remains 
anomalous, with a sign opposite to the one predicted from the single-particle density of states. The 
(N : heat capacity and non-Fermi liquid response are discussed in the context of similar results in the 

' literature. All results are obtained by evaluation of an effective single-particle free-energy operator 

D . in closed form. The method to obtain this operator is described in detail. 



o 
o 



> 



c3 



X 



I. INTRODUCTION 



• There is at present no paradigmatic 'non-Fermi liquid' in more than one dimension, in the same sense as there is the 
Luttinger liquid in the ID case. It is therefore of continuous interest to study failures of the Fermi liquid concept in 
. special cases, accessible to concrete calculation. The behavior of electrons in the presence of localized obstructions has 
long been a fertile ground for investigating deviations from simple Fermi liquid behavior. Examples are the quantum 
Hall effect pj, the Kondo problem Q, localization by impurities and the Mott metal-to-insulator transition Q. 
C< All of these save the last involve fixed impurity distributions, and their unusual behavior is related to the creation of 
I | localized states in the presence of scatterers. Because of this, states accessible to weak probes cannot be understood 
' O ■ in the picture of nearly independent quasiparticles. 

The Mott transition is unique among the above because the localized obstructions are the mobile electrons them- 
selves, due to the strong on-site repulsion between electrons of the opposite spin. This means translational invariance 
is not broken, and obstruction is dynamical, which makes the problem much more difficult than with quenched dis- 
order. Despite many contributions which have shaped our present understanding [5llr3.l^.l3.l9l ll0llllLll2j |. no picture 
of the Mott transition has emerged to date, which convincingly describes the phenomenon in space and time. 

Motivated by the Mott problem in the charge-transfer limit, which is of primary interest for high-temperature 
superconductors |13| , the present work describes non- Fermi liquid behavior in one model of the Falicov-Kimball [l4j 
group, which are intermediate between the quenched-disorder case and the fully dynamical Mott case: the disorder is 
, heavy, but annealed. In this way translational invariance is restored at the level of the ensemble, while the dynamics 
still relates to static disorder. The present model is the only one in the group in which the disorder is truly geometric, 
i.e. quantum processes do not physically interact with any classical degree of freedom. This makes it an interesting 
test-bed for non-Fermi liquid behavior, and indeed it is found that geometric obstruction induces collective behavior in 
the 'strange metal' state beyond half-filling. It somewhat resembles normal 3 He at 0.5 K, where the entropy saturates 
at In 2 per atom , simply because the atoms push against each other, obstructing kinetic motion. In the limit of 
low hopping overlap, the model develops an unexpected similarity to the Kauzmann paradox |l6j| in vitreous liquids, 
. despite the fact that all calculations are at equilibrium, so one cannot speak of kinetic slow-down in the usual sense. 
For intermediate-to-large overlaps, the model appears at first sight as a renormalized Fermi liquid. However, the 
'equilibrium thermopower,' d[ijdT, behaves exactly oppositely to what is predicted from the effective one-particle 
O ■ density of states (Tp-DOS) at the same temperature and filling. This persistence of quantum collectivity is observed 
because geometric obstructions take the form of boundary conditions, so they cannot be suppressed by hopping 
fluctuations. 

The technique used to analyze the equilibrium properties of the model is presented here for the first time. It 
enables one to find a closed analytical expression for the free energy. The essential idea is to remain in the Fock 
space of the mobile fermions, while evaluating the logarithm of the annealed partition function. The one-particle part 
of this logarithm is then itself a one-particle operator in Fock space, and can be diagonalized exactly, just like any 
tight-binding Hamiltonian. With such a careful treatment of Pauli correlations, it is sufficient to expand the annealed 
free energy to the first non-trivial term in a Kubo cluster cumulant expansion 17]. The resulting analytical expression 
is exactly correct in both the atomic and metallic limits, and qualitatively compares well with numerical calculations 
on similar models. 
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FIG. 1: The random-tiling model. Black lines are hopping integrals for up-spins, connecting copper (black) and oxygen sites 
(white). The hopping integrals are turned off where there is a down-spin on a copper site. 

II. THE MODEL 
A. Brief description 

The geometric Falicov-Kimball model used here has been introduced previously [l^. ft is given by the following 
modification of the well-known Emery (three-band) model |l9j: 

H = £d^ri i( T + e p ^2(n xM + n yAa ) + t^(l - %) (d^P^ + ^T^t) ' W 

where P ia = , 2 CT + / 2 CT — Pi+x/2 a ~ Pi+y/2 a describes the four oxygen p-orbitals around a copper d-orbital, 
and hi a , n x ,ia, n y ^ a are the number operators of the copper and oxygen sites. Notice that only up-spins can hop. 
This is a simple realization of the idea behind Gutzwiller's approximation: fermions of one spin orientation see those 
of the other 'as if occupying a band of width zero' 6]. The static down-spins are scattered at random. If one lands 
on an oxygen site, it pays the price of a charge-transfer energy A p d = s p — e& > between oxygen and copper, but is 
otherwise inactive. Thus the mobile up-spins see a lattice with defects, different for each arrangement of down-spins. 
The situation is summarized in Figure^ (For those familiar with the Emery model, note that the up- and down-spins 
here correspond to physical holes in the copper d-orbitals, not electrons.) 

In this variant of the Falicov-Kimball 01 model, there is no physical scattering off the heavy classical spins. Instead 
the mobile spins see a 'torn' lattice with open boundary conditions where the hopping has been turned off. Thus they 
are in the strong-coupling limit on the one hand, and there is no singularity in their Hamiltonian, on the other. Now 
there intervenes a significant topological simplification, namely the 'crosses' of hopping overlaps impinging on copper 
sites tile the plane, as visible in Fig.^ This gives the model its name, random-tiling (RT) model, and makes it possible 
to express the annealed free energy of the problem by a Kubo cluster-cumulant expansion without leaving the Fock 
space of up-spins 0,H(j. Double occupation of copper sites is disallowed by a separate factor exp(— (3UJ2i "iT^l) m 
the partition function, with U 3> A pc ;. This factor is not entangled with the Hamiltonian, so it is not an interaction, 
but pure bookkeeping, to force all the mobile spins away from the 'torn' part of the lattice. The projector (1 — hii) in 
Eq. JIJ is treated exactly, which is essential [2l| : both the Pauli principle and the geometric obstruction are treated on 
an equal footing. One projects the Fock-space Kubo expansion onto the single-particle subspace, and keeps the first 
non-trivial term, while higher single-particle terms can be shown to be negligible. This gives a one-particle quantum 
free-energy operator, whose spectrum can be exactly obtained in closed form. It is built from up-spin wave functions 
which anneal the geometric disorder. It provides the effective background for all residual (two-particle) interactions, 
so it is of interest to study by itself. The Mott-Hubbard transition has already been studied in some detail 22] . 

B. Detailed derivation 

The present section is purely technical, and quite independent of the rest of the article. It describes how the 
annealed free energy of the random-tiling model can be written without leaving the Fock space of up-spins. The one- 
particle part of this expression takes the form of a rapidly convergent cluster cumulant expansion around the atomic 
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limit. The main strength of the procedure is that it never violates either the geometric or the Pauli constraints. In 
particular, it is shown that the same closed expression exactly recovers the metallic free-fermion limit for the up-spins, 
when the down-spin concentration goes to zero. This is the opposite limit from the atomic one, which attests to the 
advantage of working in Fock space. Such physical considerations are collected in the last part of this section, to 
which one may prefer to skip at first reading. 

1. Problem and method 

The physical picture described above corresponds at equilibrium to the evaluation of the canonical partition function 

Z = tr e- pH (2) 

where H is the Hamiltonian fll]). under the geometric constraint, that mobile up-spins do not land on the copper sites 
occupied by the down-spins, to which the hopping has been turned off by the projector in As already mentioned 
above, this is implemented by actually evaluating the unconstrained trace of 

Z = tr e^ H e'^ v = tr e -W 2 e-^ e ^ 2 , (3) 

where the disentangled term e~@ v enforces the constraint, with 

V = Uj2^nn. (4) 

i 

(The form with V/2 turns out to be more convenient for later manipulation.) It puts the mobile up-spins onto that 
part of the lattice which is not 'torn' by the down-spins on copper sites, acting through the projector 1 — hii in 
Eq. One can easily put U — ► oo in the final expressions, but it is neater to keep a finite U 3> A p d- As expected, 
the results are numerically independent of U in this limit. 

The trace over down-spins in Eq. [21 is trivial to perform, and results in rather laborious expressions of the generic 
form 

£ f(A ll+ ... + A lM )^U LM . (5) 

{i 1 ,i 2 ,...,iAi}CC 

Here Ai are site-labelled operators in the Fock space of up-spins, essentially the crosses in Fig. This is obvious, 
since a given configuration of down-spins corresponds to some choice of M lattice sites where the hopping survives, 
out of a total L sites (L is volume) comprising the lattice C. Note that / is a function of Fock-space operators for the 
up-spins. It is something like 

f(X)~ ex p[-p(H +X)}, (6) 

where Hq is the site Hamiltonian. The Ulm are the thermal evolution operators for the annealed disorder with 
M hopping crosses present. The main point of the method used here is to reorganize expressions like (|SJ) by a 
combinatorial inversion which is valid for arbitrary non-commuting operators. This allows one to extract the one- 
particle part of the free energy expansion exactly, and show that corrections to its first (one-site) term are controlled 
in the range of fillings of interest here. 

Physically, the method boils down to a combinatorial interpretation of Kubo's cluster cumulant formula for the free 
energy. This is explained at length in Ref. |2C| . where these ideas are used to cover familiar ground. In particular, the 
reader is reminded that the expansion parameter . co rresponding to the concentration in the classical virial expansion, 
is an occupation probability in the general case [l7| . In the present work, one implements Kubo's formalism by first 
performing the combinatorial inversion at the level of the stochastic evolution operators JSJ , and taking the logarithms 
only in the last step. To be explicit, introduce the binomial transform 

u l {p) - J2^r- M ( L p Z^)u LM . (7b) 



M 
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When one of these expressions is inserted into the other, the result is an identity, so they are valid independently of 
what the Ulm may be. Then a very short exercise shows that the generating function for the C/lm's may be expressed 
in terms of the Ul(p) as 

E = (! + X ) L E V PU l{p), (8) 

M p 

where y — x/(l + x), < y < 1. This result depends only on the relation (|7a|) . and in particular, nothing is assumed 
about commutativity. Since the U lm are canonical evolution operators, the sum in JSJ is the corresponding grand 
canonical operator. Hence x involves the chemical potential for the hopping crosses, x = exp(/3i/), say. In the present 
context, it will be more convenient to express the expansion parameter y in terms of the chemical potential of the 
static down-spins, call it v. Since the hopping cross is absent when the down spin is present on a copper site, the two 
chemical potentials are effectively connected by a particle-hole transformation. When the energy origin is conveniently 
chosen so that Ed — 0, this amounts to v = —v', i.e. 



1 + e' 31 ' ' 



For the range of concentrations considered in this article, half-filling and beyond, the chemical potential of the static 
down-spins will turn out to be at the oxygen energy, v « e v > 0, so that y is exponentially small. Thus it provides 
rapid convergence in the expansion of the one-particle part of Eq. for these concentrations. 

2. Main expressions 

Set the bare copper energy Ed = 0. In Emery-model jargon, the 'hole picture' is assumed, so e p > Ed = 0, and as 
mentioned above, U 3> e p . Call 

L 

i=i 

the number operator on the copper sites, which appears in ftp. The actual forms encountered in tracing J3J are 
slightly more complicated than suggested in the previous section, because of the side factors with V/2, and because 
the down spins can also land on the oxygen sites. Overcoming these complications by ordinary diligence, the Ul(p), 
required in (JSJ, are given by 

U L (0) = e -/%^T+tw dT ) j (n) 

L 

U L (l} = ^ ^ e -/3(^/2)(^d T -«u) e -/ 3 [^^PT+t(4T P .T+ P /T d .T)] e -^^/ 2 )(* d T-™.T) „ f/ L (0)^ . (12) 



Here, N p f is the number operator on the oxygen sites. Stopping at first order in y in the general expansion (jSJ), the 
operator to be traced over the up-spins now appears as 



E 

Mi 



qL 



e -fi ep (Ni-M L ) 



[[x L - M >]] (1 + x) L exp {[In (U L (0) + yU L (l))] lp } , (13) 



where the notation '[[a; r ]] . . .' means 'coefficient of x r in . . .', and 'lp' means 'one-particle part.' Here iVj is the total 
number of static down-spins, and M± is the number of down-spins which have landed on the coppers. The number 
of oxygen sites per unit cell is q = 2. In the next section, the logarithm of Eq. (|13f) will be explicitly computed in 
Fock space. Assume for now that this has been done, and that the one-particle operator in braces in l|13f) has been 
diagonalized. Then the trace can also be performed in the space of up-spins. Introducing a chemical potential fj, for 
these latter, the trace now reads 

^\N,-MJ e A U JV Te -/3 I /(£--MJ ' y } 

Mi v 1 l/ 
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where x has been called e P v , v is the chemical potential of the down-spins, and 

ZQ*,v) =Jj(l + e /, <"- e «< fc '"») (15) 



is the grand partition function in up-spin space, with e a (k, v) the energies of the up-spin normal modes, indexed with 
a band index a, and depending on the chemical potential of the down-spins, v, as well as on a wave vector, k. Finally, 
the sum in (|14fl is evaluated by taking its largest term, resulting in the following model: find the extremum of 



,.— (l + e ^)(l + e ^-^)^^Z(^) 

V^effUi, v) = eP^eP^i ' ( ] 



in the space of /i and is, for given concentrations n-f and n^. The extremal equations are 

1 dhxZ{n,v) 



'"I 



(17a) 



^ e^-^) . 1 dlnZ(^) 

»i = TT^ + g IT^(^) + I ^ ■ (17b) 

Eq. (|17b(l is physically transparent: there are static down-spins, sitting in one level at Ed = 0, and q — 2 levels at the 
oxygen energy e p . The last term, affecting the count of down-spins, comes from the interaction with the up-spins, as 
coded by the dependence of their dispersion on v. This dispersion is derived in the next section, giving a complete, 
closed model. Its dependence on the chemical potential of the down-spins is expected on general thermodynamic 
grounds, since the presence of down-spins affects the configurations available to up-spins. 

3. Extraction of the one-particle part 

In the previous section, the binomial transform was used to treat the stochastic evolution operator induced by the 
geometric constraint (annealment). The manipulations were mostly an exercise in bookkeeping. In this section, the 
Pauli constraint is considered, and the main tool is the evaluation of analytic functions of matrix arguments. This is 
the straightforward way to evaluate analytic functions of Fock space operators, if one is only interested in the resulting 
one-particle part. Analytic functions of matrix arguments are evaluated as usual, by diagonalization: 

F(A) = FiSS^ASS- 1 ) = SF(S- 1 AS)S~ 1 . 

To prepare the ground requires, however, one last combinatorial trick. Namely, take the f(X) in JSJ to be ln(J7i(0) + 
yX). Then ln({7i(0) +j/£/l(1)) can play the role of an 1 Ull in 10), since Ul(1) in iff^) is a sum over all site-indexed 
terms, each playing the role of an 'A^ in JSJ. The expansion l|7a|l with M — L then gives, for this case, 



-/3(e p N P i + UNdi) 



+ J2 On [Ul(P) + yB l ] - In [U L (0)]} +£\ . . , (18) 

i=l i<j 

where the sum over i < j is already a second-order term. Here Bi denotes the i-th term under the sum in IJ12H . The 
effective free energy is then 

1 L 

F eS = e p N p] + UN d] -} lp , (19) 

where the braces refer to the same expression as in Eq. (|18|) . Higher order terms are dropped, because they correspond 
to higher powers in y. 

From now on the way is clear, provided one has an algebraic manipulation program. The effective free energy Ijl9(l 
may be rewritten 



F eS = to. (l + V?[U +y(U 1 -U )}V i 



i P 

-V//diag( £p /2, e p /2, 0, e p /2, e p /2)V tl (20) 



6 



where Vi = c °l(Pi-£/2 pPi+y/2 f ^'^-5/2 ^ s ^ e commn °f Fock-space operators, referring to atoms 

centered on the i-th site, with d ia the copper operator. Explicitly, Uq = diag(e~' 3 ' rp — 1, e _,3ep — 1, e - ' 3 ' 7 — 1, e _,3ep — 
l,e-^ -1), and 



U l -U Q = 



with 

TO 

4r 2 

where the following shorthand has been used: 



/ r 2 -T2 Tl r 2 -t 2 \ 

-T 2 T 2 — Tl -T 2 T 2 

Tl -Tl T Tl -Tl 

T 2 — T 2 Tl T 2 — T2 

V -T2 T 2 -Tl -T 2 T 2 



2 V , P • 2 V 

u; cos — H sin — — c, 

2 w 2 



-wsinyj - 



■ sm(p, 



• 2 V P i <P 

w sin 1 cos p, 

2 w 2 y ' 



w 



p J3A pd ^ p _ p -Pe p ^ c 



-0U 



and 



2i 



sin y> = 



COS if 



e p /2 



while the exponent in w; is 



A p£ 2 = 



(I 



■4t 2 -^. 
2 



(21) 



(22a) 
(22b) 
(22c) 

(23a) 

(23b) 
(23c) 



Coming back to F e g , the one-particle part of the logarithm may be written explicitly, giving 

L 



(23d) 



(24a) 



where 



T 



( 


£p _i_ t 

2 ^ l VP 


tpp 


tpd 


tpp 


tpp 




tpp 


^ 4- t 
2 ^ l VV 


tpd 


tpp 


tpp 




tpd 


tpd 




tpd 


tpd 


V 


tpp 


tpp 


tpd 


2 ' 6 PP 


tpp 


tpp 


tpp 


tpd 


tpp 


fp _|_ + 
2 ' PP 



(24b) 



Before the expressions for the matrix elements of T are given, note that (|24|l is like a simple tight-binding Hamiltonian, 
which may be diagonalized by a Fourier transform. The spectrum consists of three bands, one of which is nonbonding: 



'e p - As d 



£ = e P , (25b) 

with 7 2 = sin 2 k x /2 + sin 2 fc a /2. The dispersion 125|) is the effective spectrum of up-spins in the presence of static 
down-spins. It was denoted, schematically, as e Q (fc, v) in the definition l|15|) of Z(/i, v). 
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In order to specify the model (|16fl in full, it only remains to write down the matrix elements of T in i|24|l . Define 
C± by 



± 



p — c 



(26) 



where y = 1/(1 + e^ v \ and the rest of the notation is as in I|22I23|) . while v is, as always, the chemical potential of 
the static down-spins. Then 





c+ 


-C- 

4 


4£pp — 


c+ 


+ C- 

2 


Ae d = 


C+ 


+ C- 

2 



sin^, 



where 



simp 
cos tp 



2yn 



■ COSIp — £ p , 

■ cosip, 



(27a) 
(27b) 
(27c) 



[ 2 T £ +J/(2r 2 



in 

2 , 



4y 2 T 



2^2 



[^+y(2r 2 -f)]+4y^ 
This completes the derivation of the effective spectrum (|25f) of the model (|16|) . 



4- Physical considerations 



Let us first summarize the physical truncations made in the preceding derivation. Two distinct meanings of 'first 
non-trivial term' were employed. One is to keep to the one-particle part of the Fock space. This is guaranteed by the 
Fock-space formalism itself, or more technically by the fact that one has passed from operators to numbers only in the 
last step of the calculation, when diagonalizing the operator l|24|l . The other is to keep to first order in the expansion 
in y. This is a little more tricky, because the final spectrum l|25|) of the free energy is non-linear in y. Here 'first-order' 
actually means first-order at the level of the evolution operator. This is especially evident in the forms l|19|l and H20JI . 
Operationally, the rule is to keep to first order in y under the logarithm, but the one-particle part of the resulting 
logarithmic operator must be extracted exactly, i.e. without further truncation in y. 

The above one-particle effective free energy is formally the first (one-site) non-trivial term in an expansion around 
the atomic limit. In terms of the expansion parameter y = 1/(1 + e ,3ly ), where v is the chemical potential of the 
down-spins, this is an expansion around y = 0, when the copper sites are completely occupied by down-spins. It will 
now be shown that the result is correct also for y = 1, which means no down-spins at all (y — > — oo). Then the exact 
non-interacting band is recovered. This is a very useful check that the two truncation prescriptions discussed in the 
previous paragraph are mutually consistent. In particular, it shows that keeping to first order in y does not violate 
the Pauli principle. This is a benefit of working in Fock space — in similar expansions around the classical limit, the 
Pauli correlations are only taken into account in an inclusion-exclusion procedure [2Jj | , which remains approximate to 
all finite orders in the expansion parameter. 

Taking v — > — oo in Eq. (|2*S|) . one gets 

cos?/' — *• — cosy, sin ip — > — sin ip. (28) 
Putting all this together, the hopping amplitudes (|27|l become 

t P d — * t, t pp — > 0, Asd — * 0, 

which gives the dispersion for the free limit of the Hamiltonian Q. This confirms that the one-particle part of the 
original problem has been treated correctly. It is also easy to verify that putting y = gives the same result as t = 
in the original Hamiltonian, but this is of course expected. 
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FIG. 2: Entropy of mobile spins for t = 0.25 eV (left three curves), and t = 0.5 eV (right four curves). Curves are marked 
by the corresponding concentration. Saturation plateaus are clearly visible. Here A p d = 3 eV, U = 10 eV, and n = 2n-f = 2ri[. 



It will now be shown that the one-particle part of the expansion in y is controlled, for the concentrations considered 
in this article. The higher order terms in y contain one-particle contributions coming from the diagonalization of 
larger 'molecules,' such as CU2O7 (say), which appears when two 'hopping crosses' share an oxygen. This has the 
typical structure of a connected-cluster expansion. The general term is y p UL(p). For down-spin concentrations beyond 
half-filling, the saddle-point equations give v « e p , i.e. y ~ e~ /3£p = e~P Apd (remember we put Ed = 0). Note from 
Eqs. (|20[1 and 1] -?•!!) that the magnitude of the term multiplying y is roughly w = e +^ A pd j so the overall estimate for 

the first term is ylJ^l) ~ yw ~ e~P( Apd ~ Apd> < 1. Physically, A p d is just the hybridization energy associated with 
the diagonalization of a single CUO4 'molecule.' A larger cluster will have a larger hybridization gain, but it will also 
be suppressed by a larger power of y. Since hybridization gains quickly saturate as a function of cluster size, they 
cannot offset the addition of a whole bare A p d with each new Cu site added to a 'molecule.' Hence the expansion is 
rapidly convergent in the range of fillings considered below. 

For down-spin fillings below the Mott transition, i.e. corresponding to the lower Hubbard band (LHB), v w Ed — 0, 
meaning y is some appreciable fraction of unity, not exponentially small. Then presumably higher-order single-particle 
terms could play a role. This is consistent with the fact that in the LHB, there is a significant contribution of the 
interaction to the entropy, while the same is exponentially small beyond the Mott transition — the system goes from 
liquid to gas, absorbing latent heat |lg,|22|. Nevertheless, the fact that the exact band result is recovered in the limit 
y = 1 suggests that the present calculation should not be qualitatively incorrect even in the LHB, not of main interest 
here. 



III. RESULTS 



A. Small hopping overlap 

In the model, the entropy of mobile spins is formally given by the usual single-particle expression 

St = - / lU(s-)HU(s-)) + (i - //.(£-)) Mi - U(e-))] , (29) 

JBZ 

where fj, is the chemical potential of the mobile spins, while the effective bonding band e_, given by Eq. (|25[). is 
a function of v, the chemical potential of the static spins, and of the temperature. Both /i and v change with 
temperature, since the particle number is fixed for the up- and down-spins separately. In Fig. [2 the entropy of mobile 
spins is shown as a function of temperature, for two values of the bare hopping and several concentrations. The most 
striking feature of the figure are broad plateaus, or inflections, indicating dips of width ~ 100 K in the heat capacity 
cy of the system. The curve at n — 1.8 even shows a 'Kauzmann paradox' |16|. where extrapolation of the high-energy 
part would give a negative entropy at zero temperature. Of course, the entropy is strictly zero at T = 0, since the 
free energy (I16fl is microscopic by construction. The plateaus correspond to regions of low cy, associated with the 
difficulty of finding new states, which gives peaks in cy (not shown) between T = and the plateau positions. 

The plateaus are identified with static disorder in the obstructed copper sites. The proof is as follows: call the 
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T[K] 

FIG. 3: Comparison of Sf (T) (thick line) and Sf t cu{T), Eq. 1301 (thin line), for t = 0.31 eV, n — 1.2, and other parameters as 
in Fig.H 
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T[K] 

FIG. 4: Chemical potential of mobile spins (thick line) corresponding to S-\ (T) in Fig. |3] The bottom and top of the effective 
band are given by thin and dashed lines, respectively. 

concentration of up-spins on the coppers m-f . If they were static, their contribution to the entropy would be just 

St, c« = —m-\ ln(m T ) - (1 - m T ) ln(l - m T ), (30) 

and beyond half-filling, each copper site would be singly occupied, either by an up-spin or a down-spin. Thus 
TOf = 1 — m|, in the static case, where is the concentration of down-spins on the coppers. Since mj, is easily 
obtained from the calculation, one can use it to construct S-\ t cu{T) from Eq. (|30J) . and compare these values with the 
actual Sf(T) from the same calculation, given by Eq. (|29|l . As shown in Fig. |3J the hypothetical static entropy i|30|) 
coincides with the calculated entropy at the inflection point, beyond which the latter begins to rise. To show this 
more clearly, the hopping overlap is slightly increased in this figure, so the plateau is not quite flat. 

This means that in the regime t -C A pc j, excitation of the system proceeds in two steps. First the copper positional 
degrees of freedom are fully saturated: their entropy cannot be greater than that of the correspondig static disorder. 
Then the translational ones, involving the oxygens, are activated, but only after a small (~ 100 K) region is crossed, 
where there are no new states. The two steps are better separated when the bare hopping is smaller: note the 
difference between the curves for t = 0.25 eV and t = 0.5 eV in Fig.0 

Figure 0] shows how the model technically achieves this behavior. As usual with single-particle effective models, 
collectivity enters through the solution of the saddle-point equations, which fix the chemical potential at a given 
concentration and temperature. The chemical potential traverses the band fairly rapidly as the temperature increases, 
and leaves it roughly around the temperature at which the plateau appears. As it passes through the region where 
the lp-DOS drops, there appears a lack of new states for the system. Once in the (Maxwellian) tail of the Fermi 
distribution, the oxygen degrees of freedom are further excited. In this sense, the system is again analogous to 3 He. 
However, there is no regime of classical kinetic motion in the model, since the bands are built-in by assumption. 
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FIG. 5: Left column: entropy of mobile spins 12911 . Right column: equilibrium thermopower 1311 . Upper row: rif = nj, = 0.6. 
Lower row: rif = 0.6, ny = 0.001. Lines, thin to thick: hopping t = 0.25, 0.5, 0.75, 1.0 eV, respectively, with A p d = 3 eV. 



B. Large hopping overlap 

Collective effects persist in the model even when the hopping overlap becomes so large that the effective band 
appears to be normally metallic, with a well-defined Fermi surface. They are no longer immediately evident at the 
level of the entropy, but remain in the susceptibility. This is shown in Figure For comparison, the lower row shows 
the limit of free fcrmions. The entropy for = — 0.6 (upper left panel) rapidly begins to look 'normal' as the 
hopping overlap increases (thicker curves). 

The right column in the figure shows that this impression is misleading, however. It gives the equilibrium suscep- 
tibility, corresponding to the thermopower, 



In a simple rigid band, S eq should change sign from positive to negative as the filling goes beyond one-half. This 
is because the chemical potential goes down with the temperature, if the lp-DOS goes up with the filling, and vice 
versa. 

For the free model (lower right panel), the equilibrium thermopower (|31|l is indeed negative. Now add down- 
spins until the model undergoes its Mott-Hubbard transition . Then the up-spin S eq (|31f) becomes positive, even 
though the Fermi surface does not change with respect to the free situation. It is physically clear what is happening: 
the kinetic blocking effects are strongest around half-filling, so going beyond it increases the available degrees of 
freedom, since the mobile up-spins on the oxygens are not blocked. A similar outcome is expected from the effect of 
antifcrromagnetic (AF) correlations above the transition, however there it would be traced to the appearance of a dip 
in the single-particle spectrum. Namely, such a lp-DOS has a minimum around half-filling, so that beyond it, it rises 
again, just like at the bottom of the band. Thus in the spin channel, an anomalous susceptibility is due to a normal 
Fermi-liquid response. By contrast, in the charge channel and beyond the Mott-Hubbard transition, the response 
cannot similarly be inferred from the single-particle spectrum. Even though the model i|lb|) is formally single-particle 
and in the metallic regime, it is not a normal Fermi liquid. 



C. Heat capacity 

The heat capacities of the random-tiling model are shown in Fig. [fj] (thinner lines are always closer to half- filling) . 
In the LHB (left column), there is a transfer of spectral intensity with doping, between a low- and high-energy scale, 
both of which already exist at half filling. The latter is still significantly lower than the bare charge-transfer scale 
A p( j, which takes the role of U in the present context. This is an effect of hybridization, as mentioned above: the 
larger the hopping overlap, the lower this 'high' scale becomes (note that the horizontal axes in the two rows are 
different). Below the Mott-Hubbard transition, the equilibrium thermopower is of the same sign as in the free case 
ni — > (Fig. 0, even for 2n T = 2n l = 0.9. 

Concentrations beyond half-filling cannot be studied in a one-band model. In the present three-band setting, there 
appears another low-energy peak between the two peaks present at half-filling, which quickly draws strength from the 
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FIG. 6: Heat capacities of the random-tiling model. Left (right) column: below (above) half-filling. Thin to thick: left, 
n = 2n T = 2n l = 0.9,0.7,0.5,0.3,0.1, right, n = 2n T = 2n L = 1.0,1.2,1.4,1.6,1.8, respectively. Upper row: t = 0.5 eV, 
A p d = 3 eV. Lower row: t — 1 eV, A p a — 3 eV. Vertical dashed line: renormalized charge-transfer scale A p< j, Eq. l23dl . Insets: 
upper, detail of the narrow peak; lower, n = 1.4, t = 2 eV, A p d = 3 eV. Black triangles on the axes: W/t, Eg. 1341 




high-energy peak as the doping increases. The new low-energy scale induced by doping is specific to the three-band 
model in the 'in-gap' concentration range. All the relevant scales can be extracted analytically in the limit T — > 0, 
since in that limit the chemical potential of the down spins is forced to one of two known values: v — > e d = in 
the lower Hubbard band (LGB), and v — > e p in the in-gap band (IGB), beyond the Mott-Hubbard transition. The 
expressions (|27[1 are easy to evaluate in these limits. Below the transition, the width of the lower Hubbard band is 
practically unchanged from the non-interacting value 

"2 A 

(32) 

where A pd = e p — e d . Here the only effect of strong correlations on the band is a small but finite shift of the copper 
level downwards, 

hm Ae d = A pd sin 2 | - A pd cos 2 | < (LHB) . (33) 

The IGB is, on the other hand, strongly affected by correlations. The copper level is shifted upwards, and the band 
is significantly narrowed. Its width in the zero temperature limit is |22| 

limW= ^SL^ — (IGB). (34) 

t^o i + A pd /(2A pd ) 

The copper level shift is now large and positive: 

hm Ae d = A pd -A pd cos 2 |>0 (IGB). (35) 

Instead of using this formula, it is much more convenient to note that the distance from the bare oxygen position e p 
to the middle of the IGB is just A pd in the same limit, i.e. the uppper edge of the IGB is given by any of the two 
formulas 

lim(e d + Ae d )=Km(s p -A pd + ^\ (IGB). (36) 

The relevance of these scales for the various peaks is analyzed in Fig. [7\ Both below and above the Mott-Hubbard 
transition, the lowest-energy peak is due to intraband transitions. The peak scales perfectly with W for the IGB, 
and with Wq for the LHB. The high peak in the LHB is due to interband processes, namely excitations into the 
empty, non-dispersive non-bonding band at e p , which is a distance A pd + |Ae<i| + Wq/2 from the Fermi level when the 
occupied band is nearly half-filled from below. The middle (new) peak in the IGB is also due to interband transitions. 
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FIG. 7: Scales of the peaks in heat capacities for t = 1 eV. Top (bottom) row: n=0.9 (n=l.l). Lines, thin to thick: A p d = 
2, 3, 4, 5, 6, 7, 8 eV, respectively. The drawing in the middle of the top row schematically shows the scales as barred thick lines, 
left for the lower Hubbard band (LHB), right for the in-gap band (IGB). The band edges are given by thin horizontal lines, 
and the occupied levels indicated by shaded squares. 



As noted above, the middle of this band is a distance A p d from the oxygen level, giving the scale of the middle peak 
when the band is nearly half-filled from above. 

Finally, the upper (third) peak in the IGB is again due to interband transitions, but of a different kind. They 
are governed by the bare scale A p£ ;, which plays the role analogous to U. However, there is no effective band in the 
model at Ed beyond the Mott-Hubbard transition (the position Ed is denoted by a broken line for this reason). The 
high peak is probing the 'undressing' scale, at which the system effectively reverts to the atomic limit. Note that the 
two thinnest lines, corresponding to A p d/t = 2 and 3, respectively, do not scale so well by A p d in the right column, 
and simultaneously have not yet developed a third peak in the middle column. The separation of scales between the 
middle and high peak should be regarded as a sign of a well-developed charge-transfer regime A p d > t and n > 1, 
while for A p d ~ t, and n ~ 1 from above, the two scales are mixed by hopping fluctuations. The upper peak for these 
two thin lines has a higher scale than A p d, because the copper level they see is still effectively shifted by hopping 
correlations. This can be corroborated by an amusing 'quick and dirty' check (not shown), where one tries to scale 
the upper peak in the IGB by A p d + |Ae<i|, but using the value | Ae^l for the LHB. Even though this is formally 
inconsistent, it turns out that the scaling in the lower right panel of Fig. is improved, since the level shift |Ae<j| in 
the LHB is bigger for stronger hopping fluctuations (smaller A p d/t). 

The same two IGB regimes, hopping-fluctuation and charge-transfer, may be discerned in the upper and lower 
panels in the right column of Fig- El drawn to absolute scale t. In the upper right panel, the input t/A p d <C 1, and 
the middle (new) peak is higher than W and lower than A p d- Then A p d appears as a crossover scale between the 
high-energy peak and the two lower-energy ones. This is the charge-transfer regime, mentioned above. 

The hopping-fluctuation regime prevails when the new peak is at low energy with respect to both W and A p d- The 
lower right panel in Fig. [S] shows the case A p d = t, when some strength still remains in the high-energy peak. The 
lowest peak has already broadened, leaving no trace of the Brinkman-Rice effect (the scale of the lowest peak has 
grown to about 1000 K, which is the whole range shown in Fig. EJ) . The inset shows the hopping-fluctuation regime 
fully developed, where the upper peak has completely disappeared, so there are again only two peaks like in the LHB 
(left column). Both are low-energy phenomena (A p d/t ks 1.39 for the inset). For all of the right column in Fig. El 
the equilibrium thermopower stays of the 'wrong' sign, like in the upper row of Fig. no matter how large t is. The 
fact that the metal is both strange and has an ordinary Fermi surface is consistent with the fact that the unusual 
two-peaked feature in the heat capacity is observed on a scale smaller than the band-width. 



13 




D. Comparison with other work 



To compare the present approach meaningfully with others, some general remarks are in order. The approach 
described here works at present only for the equilibrium properties, a limitation not shared by some other available 
methods. It is in the thermodynamic limit, like other theoretical schemes, and unlike purely numerical ones. It is 
unique by being fully analytic. The spectrum 1251) is given in closed form, so the only numerical part of the work is 
to fix the chemical potentials /i and v by Eqs. I|17|) . like in a non- interacting band problem. 

One should also have some idea, what such a comparison hopes to achieve. The present paper regards models 
of Falicov-Kimball type as less interesting for their own sake, than as stand-ins for problems one would really like 
to treat, the dynamical Mott problem in particular. In this context there is no a priori advantage to one model 
over another, and it is preferable to concentrate on their similarities. Differences are significant only if they can be 
settled by a reliable external criterion, such as a known exact result. Hence the purpose of this section is to identify 
some 'generic' behavior of this kind of models, see how the present work fits in, and hopefully understand better the 
relationship to the Mott type of problem. The papers used for comparison here are selected as representative of recent 
work, without implication as to merit or priority relative to others. 

El Shawish et al. |24| have studied a two-band spinless Falicov-Kimball model, which may be mapped onto a 
Hubbard model in which the two kinds of spins, up and down, hop with different overlaps, tf ^ t^. The method was 
purely numerical, Lanczos algorithm on a lattice of 4 x 4 sites. The limit if = corresponds to the original model [l4j |. 
and their result for this case is shown in the bottom panel on the left in Fig. EI The full curve for half-filling is directly 
comparable with the n = 1.0 curves in Fig. [21 Obviously the general trends are the same, a plateau at S = In 2, 
dropping steeply to S — at T = 0, and rising very slowly at high temperatures, due to incoherent hopping. However, 
in Fig. the entropy curves show a gap in the spectrum, while in Fig.[5]there is only a tendency for the effective mass 
to diverge, or band-width to collapse, like in the Brinkmann-Rice transition, broadened here to a crossover. Recall 
that the first peak in the heat capacity is scaled by the band- width W , which is very small when t <C A pc ;, meaning 
the entropy plateau is pulled in to very low temperature instead of being gapped. The gap is missing because only 
the one-particle part in Fock space is kept; one expects that the two-particle terms would indeed produce a gap near 
half-filling, once the band-width in the one-particle sector became small enough. (It is interesting that the down-spin 
entropy in the present model goes to zero at the Mott tranisition [li|. indicating an ordered pattern, even in the 
absence of two-particle terms.) 

The other panels in the left of Fig. [H] show the evolution in tj/tj., up to the full Hubbard model ff = t± in the top 
panel. The evolution from small t/A pc [ to large t/A p( j in the upper left panel of Fig. [S] is remarkably similar. This 
could have been expected, because the two spin orientations correspond to the / and d bands in the original two-band 
spinless Falicov-Kimball model. Hence an equality between the / and d bandwidths is naturally similar to a small (in 
units of t) splitting of the copper and oxygen energies in the present work. The similarity nevertheless raises some 
issues in relation to the Mott problem, which will be touched upon in the Discussion below. 

Hettler et al. |25| have studied the two-dimensional one-band Falicov-Kimball model by the dynamical cluster 
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FIG. 9: Heat capacities of the Hubbard model, from Lanczos |2^| (left) and quantum Monte Carlo [27) (right) calculations, 
discussed in the text. Note that positions of panels with the same concentrations are shifted. 



approximation (DC A). This is an extension of dynamical mean- field theory (DMFT [lj) to clusters of more than 
one site, embedded in the same, spatially unresolved host medium. In infinite dimensions, the simplest, one-site 
DCA reduces by construction to the DMFT 0,[2f|. In finite dimensions, it reduces to the latter's finite-dimensional 
analogue, the dynamical mean-field approximation (DMFA). The DCA results |25( for the entropy are shown in the 
right panel of Fig.[Sl for a one-site cluster (DMFA, full line), and four-site cluster (broken line). The DMFA retains 
the full In 2 per site entropy in the ground state. The four-site DCA does better, but still appears to saturate at 
about one- half that value near T = 0. Hence the DCA agrees with the present work, and with the numerical Lanczos 
results above, in the high-temperature regime, but does not cross over correctly to low temperature. 

The heat capacities for the IGB in the right column of Fig. Uncross at two fairly well defined points. Vollhardt [2^ has 
noted this to be a quite general feature of strongly correlated systems, for families of heat capacity curves parametrized 
by pressure or concentration in experiment, or the on-site repulsion U in the one-band Hubbard model. The crossing 
of concentration-parametrized families of curves in the IGB, and their failure to cross in the LHB, fits well with 
previous observations based on Fig. [SJ that it is the IGB, rather than the LHB, which is singularly affected by strong 
correlations. 

A two-peak structure in cy is ubiquitously found in numerical studies of Hubbard-like systems 0, [25[ [2(| [27} ■ In 
Fig. [5] on the left arc the results of a Lanczos calculation [2^ for the one-band Hubbard model on a 4 x 4 lattice. The 
right-hand panel is from a quantum Monte Carlo (QMC) calculation on a 6 x 6 lattice for the same model. They 
do not cross at any one point, and have two peaks, both of which are generically features of the LHB in the present 
model. This is not surprising, since concentrations below half-filling are the only ones relevant in the one-band model, 
because of particle-hole symmetry. 

It has been noted j2|J that the two calculations in Fig. differ significantly in their outcomes. First, in the QMC 
results the first peak loses strength much faster with decreasing doping. Second, in the QMC results the peak also 
moves, while in the Lanczos calculations it stays much in the same place. From the RT model's point of view, the first 
peak is expected to scale with the non-interacting bandwidth Wq, so its position should not depend on concentration, 
except insofar as due to the changing position of the Fermi level. This is visible in the left column of Fig. [21 The first 
peak there barely moves, but it moves a little more in the bottom left panel, where the non-interacting band is wider, 
so the Fermi level shifts more rapidly with concentration. Also, the narrower the LHB, the sharper the first peak, and 
the more persistent as the concentration decreases. This suggests a physical interpretation of the differences in Fig. [5] 
The QMC heat capacities imply that the effective band seen by one-band Hubbard electrons below half-filling rapidly 
widens as the concentration decreases. The Lanczos calculation implies a situation similar to the present model in the 
LHB: band- width independent of concentration, and the effective band somewhat narrower than the QMC prediction. 

To summarize, below half-filling the present analytical method gives similar qualitative behavior of thermodynamic 
quantities as other, numerically much more involved treatments of one-band models. It is also quantitatively correct in 
the known exact limits, of high and low temperature, and high and low doping. Where it can be checked quantitatively 
in the intermediate regime, namely in Fig. [3] it is correct as well. It misses the gap in the up-spin spectrum around 
half-filling, replacing it with a BR bandwidth collapse, presumably because it retains only one-particle terms in the 
Fock-space free energy operator. The doping dependence above half-filling is qualitatively different than expected 
from the one-band analogy. At fixed Hamiltonian parameters, there is a transfer of spectral weight to a new scale, 
which emerges only with doping. The new scale is connected with interband excitations across the renormalized 
charge-transfer gap. The present calculation also parallels Lanczos calculations on the Hubbard model, both in the 
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shape and behavior of the peaks in the heat capacity below half-filling. 

IV. DISCUSSION 

The main modelling interest of the present article was to determine when the one-particle construction, adopted to 
calculate the free energy, indeed corresponds to a Fermi liquid, and when it does not. For that purpose, a minimal 
negative criterion to observe a non-Fermi liquid at a finite temperature was adopted. A system is not a Fermi liquid 
if the susceptibilities do not follow from the lp-DOS at that temperature. Put conversely, a necessary condition for 
a system to be a Fermi liquid is that the lp-DOS responds 'rigidly' to an infinitesimal change in temperature. 

It has been found previously |22j that the Mott transition in the RT model is first order, from a liquid of light 
particles, in the lower Hubbard band, to a gas of heavy ones, in the renormalized in-gap band. (While the down-spin 
entropy drops to zero at the transition, indicating AF order, the chemical potential is a monotonous function of the 
filling, i.e. there is no phase separation.) In the present work, it is shown that while the liquid is a Fermi liquid, the 
gas is not. The absence of residual interactions (which makes it a gas) has been bought at the price of an anomalous 
variation of the effective mass with the temperature, such that the thermopower susceptibility i|31|) cannot be predicted 
from the lp-DOS at any given temperature. This remains so even if the hopping overlap is large, because the model 
is always in the strong-coupling limit: changing the hopping does not affect the geometric obstructions, as obvious 
from Eq. The collectivity involved is very simple. The fermions need to agree which are the best positions for 
the annealed impurities, in precise analogy to 3 He atoms pushing against each other. 

At low hopping overlap, this geometric effect so dominates the mean properties of the gas, that it is unable to find 
new states of motion in a range of temperatures, giving rise to plateaus in the entropy. (Extrapolating experimental 
trends 0], one could imagine similar plateaus appearing in 3 He around 0.5 K, if the pressure were sufficiently 
increased.) These are remarkably similar to Kauzmann plateaus in vitreous liquids, where they are due to the 
configurational rearrangements falling out of equilibrium. Clearly the Kauzmann phenomenon |l6( is a non-specific 
sign that kinetic exploration of configuration space has been obstructed, whatever the responsible mechanism. Here 
it has an interesting connection with the Brinkmann-Rice (BR) bandwidth collapse 0, which is signalled in Fig. [2 
by a steep rise of the entropy from zero temperature at half-filling, indicating an infinite effective mass. By the same 
token, the entropy plateau means that the effective mass is zero, and it is obtained by continuous deformation of the 
BR entropy curves with doping. It may seem surprising that obstruction of kinetic motion should be associated with 
zero thermodynamic effective mass, but note that in the absence of new states, there can be neither conduction nor 
dissipation. Zero effective mass is just the equilibrium version of what is observed in the Hall effect: the longitudinal 
Hall conductance and resistance simultaneously drop to (almost) zero at the steps in the transverse resistance, precisely 
because most of the carriers are localized . 

The place of the present calculation among various cluster approaches is simple to state. It is the ordinary Kubo 
cluster cumulant expansion of the annealed free energy, with the single proviso that it is calculated in Fock space. 
The transition from operators to numbers is made at the very end, after the single-particle Fock term has been 
extracted. It is theoretically quite revealing that the results presented here could be obtained by diagonalizing a 
single Cu04 molecule, i.e. the lowest non-trivial term in the underlying cluster expansion. Clearly, the main work 
was in fact done by the algebraic and combinatorial machinery which respectively enforced both the Pauli principle 
and the geometric constraint. Once these kinematical issues were properly separated from the dynamical ones, it 
turned out the problem was dominated by the former. This outcome is physically plausible by adiabatic arguments: 
the low-energy physics of the Falicov-Kimball model is not controlled by any elaborate dynamical correlations, such 
as one would expect in the Mott case, where the scatterers and the scattered are equally light. It is also consistent 
with the previous observation that the system is a gas beyond half-filling. In the lower Hubbard band, where it 
is a liquid, higher-order clusters should affect the results, but that is not expected to be qualitatively significant. 

The comparison with other, numerically much heavier calculations in the preceding section gives some grounds 
for reflection on the Mott problem. This has been dominated for the past forty years by Gutzwiller's intuitions [||, 
that electrons of one spin see those of the other as a 'smeared background' (Gutzwiller ansatz), and further 'as if 
occupying a band of width zero' (Gutzwiller approximation). This point of view clearly implies there should be some 
similarity between Falicov-Kimball and Hubbard-model low-energy behavior. Such a similarity appears to have been 
found in one of the calculations cited above p^. in the smooth evolution of the entropy as tf jty was raised from zero 
to one. If one knew that Gutzwiller's proposals were correct, the very similar evolution of the entropy, observed in the 
present calculation as tj A pc i was increased, would be naturally expected. As things stand, one should keep in mind 
the alternative resonating- valence-bond (RVB) scenario of Anderson [l3| . in which quantum coherence between the 
two spin orientations is an essential ingredient of Mott's insulating ground state. Thus the possibility remains open, 
that the Falicov-Kimball limit acts as a kind of trap for all approximate treatments of the Hubbard model which 
incidentally destroy this coherence, over time or space. 
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It is also of some interest that the heat capacities in the RT model are qualitatively more similar to numerical 
Hubbard model results, than to one-band Falicov-Kimball results. It is possible that the oxygen degree of freedom 
plays a role analogous to the 'other' spin in the one-band Hubbard model. 

To conclude, geometric obstruction in thermodynamic equilibrium has been explored in an analytic single-particle 
quantum model of Falicov-Kimball type. In the regime t <C A pc 2, it gives rise to a plateau in the entropy of mobile 
spins. For t ~ A p( z, the mobile spins appear to be normally metallic, only with renormalized parameters. Beyond 
half- filling, their equilibrium thermopower susceptibility nevertheless reflects the microscopic kinematic restrictions, 
behaving oppositely to what is expected from the single-particle properties. For these fillings, correlation effects in 
the model do not simplify in the limit of low temperature and large hopping overlap. 
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